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ABSTRACT 

9^ Ori E is a young, moderate mass binary system, a rarely observed case of 
spectral-type G-giants of about 3 Solar masses, which are still collapsing towards 
the main sequence, where they presumably become X-ray faint. We have obtained 
high resolution X-ray spectra with Chandra and find that the system is very 
active and similar to coronal sources, having emission typical of magnetically 
confined plasma: a broad temperature distribution with a hot component and 
significant high energy continuum; narrow emission lines from H- and He-like 
ions, as well as a range of Fe ions, and relative luminosity, L^/L-^oi = 10~^, 
at the saturation limit. Density, while poorly constrained, is consistent with 
the low density limits, our upper limits being rig < lO^^cm"^ for Mg xi and 
rig < lO^^cm"^ for Ne ix. Coronal elemental abundances are sub-Solar, with 
Ne being the highest at about 0.4 times Solar. We find a possible trend in 
Trapezium hot plasmas towards low relative abundances of Fe, O, and Ne, which 
is hard to explain in terms of the dust depletion scenarios of low-mass young stars. 
Variability was unusually low during our observations relative to other coronally 
active stars. Qualitatively, the emission is similar to post main-sequence G- 
stars. Coronal structures could be compact, or comparable to the dimensions 
of the stellar radii. From comparison to X-ray emission from similar mass stars 
at various evolutionary epochs, we conclude that the X-rays in 6^ Ori E are 
generated by a convective dynamo, present during contraction, but which will 
vanish during the main-sequence epoch, and possibly to be resurrected during 
post main-sequence evolution. 

Subject headings: stars: individual (tetOl Ori E); stars: coronae; stars: pre-main 
sequence; X-rays: stars 
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1. Introduction 



Intermediate mass pre-main sequence (PMS) stars, like their massive cousins, are 
difficult to study because of their rapid evolutionary time scales. Though not as short as 
stars above 8 M0 which take less than 10^ years to reach the main sequence, intermediate 
mass stars between 2M0 and 8Mq may only take 10-20 Myr. In both cases the accretion 
time scales dominate the evolution time to the zero-age main sequence (ZAMS), in contrast 
to the low-mass T Tauri stars for which PMS contraction times are longest. Intermediate 
mass PMS stars are also not easily found and id entified, and most existing studies focus on 



Herbig Ae and Be (HAeBe) stars. Herbig stars (jHerbig]|l960l ) are recognized as such once 
they already contracted to high enough photospheric temperatures to be optically identified 
as A and B stars and are thus already close to the ZAMS. Herbig stars mark the tran sition 



between formation mechanisms of low-mass and high-mass stars (iBaines et al. 



Herb i g Ae stars seem more similar to the low-mass T Tauri stars ( iWaters &: Waelkens 



1998: 



20061 



Vink et al. 



( jPrew et al. 



2OO5I ). Herbig Be stars are more similar to embedded young massive stars 



19971 ). Both the Ae and Be stars are all already in fairly late PMS stages. 



Most studies of Herbig stars use infrared and optical wavelengths to probe their 



circumstellar disks and dusty environments. Recent studies suggest that specifica . 



HAe stars there is evidence not only for circumstellar disks (IMannings fc Sargent 



Grady et al 



ly in 



1997 



I999I ) but also indications of dust shadowing and settling indicative of dust 



grain growth and planet esimal formation ( lAcke fc Waelkens 



2004 



Grady et al.l 



pianet e 
] 20051). 



2004 



DuUemond fc Dominik 



Some studies also su ggest that magnetospheric accretio n 



analogous to classical T Tauri stars is possible (jMuzeroUe et al. 



2004 



Guimaraes et al. 



Grady et al. 



2004 



20061 ). Recent modeling of 37 Herbig Ae/Fe stars using UV spectra 



revealed that all but one show indica tions of accretion wit h accretion rates in many cases 



substantially exceeding 10 ^ M0yr ^ ( iBlondel &: Djie 



20061). 



-4- 



Binarity also seems to be an important attribute i n the formation and evolution of 
intermediate mass stars. In a sample of 28 HAeBe stars, 



Baines et al 



(120061 ) find a binarity 



fraction of almost 70% with a higher binary frequency in HBe star s than in HAe star s. HAe 



stars with close companions also seem to lack circumstellar disks (IGrady et al. 



200M). 



X-ray studies of young intermediate mass stars are still quite rare and to date 



are moderately bright in X-rays ( 


Damiani et al. 


199 


4; 


Zinnecker & Preibisch 


1994 




" 1 1 

Hamaguchi. Yamauchi & Kovama i2005l: iStelzer et al. 


2006b 


). This fact is already quite 



remarkable since main sequence A-stars lack strong winds or coronae and it suggests 
that the physical characteristics of HAe stars stars differ from those of main sequence 
A- and B-stars. Mechanisms suggested range from active accretion to coronal activity to 
some other form of plasma confinement. It is also possible, given the high frequency of 
binaries amoiig HAe Be stars, that some X-ray sources could be due to late-type companions 



( IStelzer et al 



2006a 



b|). A detailed summary can be found in a recent Chandra high 



resolution spectro scopic study of the HAe star HD 104237 in the e Chamaeleontis Group 



(^Testa et al. 



20081). 



In this paper we focus on X-ray emission from 6^ Ori E, which was recently determined 
to be an intermediate mass binary star. The Orion Trapezium is gene r ally known for its 



2001 



20031 : 



Stelzer et al. 



ensem ble of the nearest and youngest massive stars (iSchulz et al. 
20051 ). Recent studies now suggest the presence of several intermediate mass stars. 9"^ Ori 
A harbors the second most massive 0-star of the Tra pezium, but also two unidentified 



intermediate mass stars both between SMp^ and 7M0 (IPreibisch et al. 



19991 ). The system 



is particularly interesting in X-rays for its high luminosity a nd hard spectral p roperties 



as well as giant hard X-ray outbursts (iFeigelson et 



temperatures during these outbursts exceed lO^K ( ISchulz et al. 



al 



2002 



Schulz et al 



20061 ). Plasma 



20061 ). While the latter 
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authors suggested a possible link of these outbursts to binary interactions involving the 
closer intermediate mass companion, there is also some evidence that these may be 
connected to the more distant companion (M. Gagne, private communication). 

9^ Ori E is another system now known t o contain young intermediate mass stars. It was 



l ong m isidentified as B5 to B8 spectral type (jParenago 



1954 



Herbig 



1960). 



Herbig fc Griffin 



(120061 ) obtained optical spectroscopic radial velocity measurements and identified the 
system as a binary containing two G III type stars of masses of about 3-4 in a 9.9 
day orbit. Evolutionary tracks constrain the age of the system to 0.5-1.0 Myr making the 
components of 6^ Ori E some of the youngest intermediate mass PMS stars known and far 
younger than stars in the HAeBe phase. 6^ Ori E is not among the optically brightest stars 
in the Tra pezium, but has long been recognized as the second brighte s t Trapezium sourc e 



in X-rays (iKu. Righini-Cohen fc Simon 



1982 



Gagne fc Caillault 



1994 



Schulz et al 



200ll ). 



The Chandra Orion Ultradeep Project (COUP) observed 6^ Ori E (COUP 732) for a total 
exposure of about 10 days over a tir ne period of 3 week s and found a low level of variability 



including one moderate X-ray fiare (IStelzer et al. 



20051 ). Its luminosity during COUP was 



determined to be log L^. [ergs s ^] = 32.4; early Chandra High Ene rgy Transmission G rating 
(HETG) spectra indicated plasma temperatures of up to 50 MK (jSchulz et al.l 120031 ). 



The HETG Orion Legacy Project has now accumulated almost 4 days of total exposure 
of 9^ Ori E allowing for an in depth study of its X-ray spectral properties. The following 
analysis of 9^ Ori E's high resolution X-ray spectrum is aimed to characterize its coronal 
nature. The existence of coronal X-rays confirms pr edictions that yery yo ung intermediate 



19931 ) and may possess 



mass stars of less than 4M0 are not fully radiative (iPalla fc Stahlerl 
some form of magnetic dynamo. We also compare these properties with the ones observed 
in 6*^ Ori A, various T Tauri stars including the relatively massive T Tauri star, SU Aur 
(2M0), active coronal sources, and post-main sequence evolved G-type giants. The optical 
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and binary system parameters of 9^ Ori E can be found in lHerbig &: Griffin! (120061 ). 



2. Observations and Analysis 



2.1. Observations, Data Processing 



As part of the HETGS Orion Legacy Project, we have observed 9^ Ori E on 11 
separate occasions from 1999 through 2007, mostly within our HETG Guaranteed Time 
program, with individ ual exposure times ranging from about 10 to 50 ks. The HETGS 



(jCanizares et al. 



20051 ) is an objective transmission grating spectrometer with two channels 



optimized for high and medium energies (HEG and MEG, respectively). The HEG and 
MEG spectra of each point source in the field form a shallow "x" centered on the zeroth 
order image. Since the Orion Trapezium field is crowded, we had to take special care to 
avoid source confusion when possible, and to assess contamination and reject spectra when 
not. The range in spacecraft roll angles, the redundancy provided by multiple gratings and 
orders, the narrow point-spread-function, and the efficiency of order-sorting with the CCD 
energy resolution all help to provide a reliable spectrum. 



We processed the data with CIAO 3.4 (IFruscione et al.ll2006l ) taking care to fine-tune 
the zero order detection to accurately center on 9^ Ori E. Response files were made with the 
most recent cali bration database availabl e at the time (version 3.4). Further analysis was 



done using ISIS (IHouck fc Denicola 



20001), an Interactive Spectral Interpretation System for 



high resolution X-ray spectroscopy, developed especially for scriptable, extensible analysis 
of Chandra high resolution spectra. 

We give an observing log in Table (H along with ancillary information and some derived 
properties of each observation. 
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Table 1. Observation Log 





Date 


lime 


texp 


j,b 
(f) 


c 


Rate X 10 


liux xlO 


Flux'i 


xl0^2 


V 








n 1 
[ks] 






[counts s J 


[phot cm s J 


[ergs cm ^] 


n —1 1 
[kms J 


3 


1999-10-31 


05:47:21 


49 


0.55 




105 


1 105 (0 0331 


3.30 


(0.10) 


-11 (22) 


4 


1999-11-24 


05-37-54 


31 


0.93 


H 


83 


1 250 (Q 025) 


3.55 


(0.07) 


-6 (52) 


2567 


2001-12-28 


12-25-56 


46 


0.30 




69 


873 (0 0331 


2.27 


(0.09) 


65 l'55'l 


2568 


2002-02-19 


20:29:42 


46 


0.70 


HM 












7407 


2006-12-03 


19:07:48 


25 


0.33 




88 


1.085 (0.052) 


3.17 


(0.15) 


68 (49) 


7410 


2006-12-06 


12:11:37 


13 


0.60 




77 


1.030 (0.072) 


2.87 


(0.20) 


-4 (72) 


7408 


2006-12-19 


14:17:30 


25 


0.93 




78 


0.994 (0.049) 


2.88 


(0.14) 


-35 (60) 


7409 


2006-12-23 


00:47:40 


27 


0.28 




90 


1.140 (0.050) 


3.43 


(0.15) 


96 (40) 


8897 


2007-11-15 


10:03:16 


24 


0.35 




84 


1.065 (0.052) 


3.16 


(0.15) 


72 (41) 


8896 


2007-11-30 


21:58:33 


23 


0.93 


M 


25 


1.030 (0.043) 


3.00 


(0.13) 


-7 (82) 


8895 


2007-12-07 


03:14:07 


25 


0.55 




86 


1.045 (0.049) 


3.10 


(0.15) 


86 (37) 



is the Chandra observation identifier number. 

^Orbital phase {cf)) was computed from the ephemeris of iHerbig fc GriffinI (|2006l ). 

"^The "C" column indicates spectra with severe confusion, "H" and "M" respectively indicate whether the 
HEG or MEG spectrum could not be used. 

The model-independent average flux in the 2-1 7A range, computed from HEG and MEG counts, flux- 
corrected using the responses. 



'The line-of-sight velocity, heliocentric correction applied. Values in parentheses are the 90% errorbar 
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(~ 1.6a). See Section [2.4.11 for explanation. 
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2.2. Source Confusion 



To assess confusion in detail, we used two techniques. For field sour ce zeroth-order 



coinc idence with the diffracted spectra of 6^ Ori E, we used the COUP fiGetman et al. 



20051 ) source list, and for each observation we transformed the list's celestial coordinates to 
the diffraction coordinate system for 6^ Ori E. We did not find any source on the spectral 
regions with significant counts. 

The second technique assessed the contamination from sources near the zeroth order, 
such that their diffracted spectra would overlap with the HEG or MEG spectra of 6^ Ori E, 
and so close to the zeroth order that the order-sorting by CCD energy would distinguish 
orders. For this, it was crucial to inspect the events' distribution as selected from the 
Ori E default binning region (cross-dispersion region half width of 6.6 x lO^^deg) in 
diffraction distance versus energy (the CCD blurred energy) coordinates in which zeroth 
orders appear as vertical distributions and diffracted photons as hyperbolas. Here we found 
significant contamination for a few observations and had to reject all or some orders. 

The useful exposure from which we can extract spectra, light curves, and line fluxes 
totals to 260 ks. Rejected orders or observations are flagged in Table [U A cumulative 
counts spectrum is shown in Figure [TJ 



2.3. Light Curves, Variability 

There was little variability of any significance within any observation. The MEG rate 
was about 0.05 counts/s. Variations within each observation were consistent with statistical 
uncertainties — no abrupt increases or slow decays characteristic of coronal flares occurred. 
Observation to observation, there was also no variability, except for one which had a 
significantly lower count rate than the others. We show the mean flux rate per observation 
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Fig. 1. — Here is an overview of the cumulative spectrum (about 260 ks) for combined MEG 
and HEG flux, for a bin size of 0.02 A. Some significant lines are marked. The statistical 
counting uncertainty is shown in gray. The lower panel shows the residuals against the 
emission measure model. The integrated flux (1-40 A) at Earth is 3.3 x 10~^^ ergs cm~^s~^ 
(1.3 X 10"^photcm~2g-i)_ 



for MEG and H 
the ephemeris of 



G spectra from 2-17 A in the upper left graph in Figure [2], phased using 



Herbig &: GrifEnI (120061 ) ; values are also listed in Table [TJ 



2.4. Spectral Analysis 



Spectral analysis is an iterative process. For detailed spectral diagnostics and 
models, accurate line fluxes, centroids, and widths are fundamental to determinations of 
temperature structure, abundances, and dynamics. However even at HETGS resolution. 
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we are dependent upon plasma models for accurate estimation of the continuum and for 
assessment of blen ding. We base our spectral models on the Astrophysical P l asma Emission 



;d: 



Smith et al. 



2001). ionization balance of 



Mazzotta et al. 



fll998h . and 



Database (AP 

abundances of I Anders fc Gre vessel (jl989l ). We begin by fitting a one-temperature component 
model to the short wavelength continuum, then add one or two more components to get 
a reasonable match to the continuum throughout the spectrum, ignoring strong lines in 
the process. Next we fit about 100 lines parametrically with unresolved Gaussians folded 
through the instrumental response, using the line-free plasma model for the local continuum. 
To improve the statistic per bin, we group the spectra by 2-4 bins and we combined the 
MEG first orders and combined the HEG first orders, both over all observations, then 
fit the MEG and HEG jointly. (The default binning oversamples the HEG and MEG 
resolutions by a factor of two, or 0.0025 A and 0.005 A, respectively.) Combination of 
spectra is done dynamically — each effective area and redistribution matrix are distinct, 
and summed counts are compared with the summed folded models to compute the statistic. 
We adopted an int erstellar absorption column of Nh = 2 x 10^^ cm~^ as determined by 



Schulz et al. 



(I2OOII ). For some weak lines, we froze the wavelength at the theoretical value. 



This allows us to obtain a limit on the flux which can provide important constraints on 
emission measure and abundance reconstruction. 

After we have mean line fluxes and centroids for a variety of elements and ions (see 
Table [2]), we can reconstruct the emission measure distribution assuming that the plasma 
has uniform abundances and is in collisional ionization equilibrium. We use a uniform 
logarithmic temperature grid and minimize the line flux residuals by adjusting the weights 
in each temperature bin as well as the elemental abundances. Since this is an ill-conditioned 
problem, we impose a smoothness on the emission measure by using its sum-squared second 
derivative in a penalty function. 
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Emission measure reconstruction is also an iterative process. First we ignore lines with 
large wavelength residuals relative to their preliminary identification based on expectations 
of a baseline plasma model, since they are likely misidentified or blended. Then we 
reconstruct a trial emission measure distribution. Lines with large flux residuals are 
rejected, since they may be symptomatic of unresolved blends or inaccurate emissivities 
due to uncertainties in the underlying atomic data. We repeat the fit with the accepted 
lines. We use the emission measure and abundance model to generate a synthetic spectrum 
and compare to the observed spectrum. Here we can adjust the line-to-continuum ratio 
by adjusting the normalizations of the emission measure and relative abundances. If the 
continuum model was improved, we start over by fitting the lines with the improved 
continuum. Finally, we perform a Monte-Carlo series of fits in which we let the measured 
line fiux vary randomly according to its measured uncertainty. This provides an estimate of 
the uncertainty on the emission measure and the abundances. 



- 13 - 




Fig. 2. — Left panel, top: the circles show the mean fluxes of the HEG and MEG spectra over 
the wavelength range 2-17 A for each unconfused observation (see Tabled]). Error-bars show 
3(7 uncertainties. Left panel, bottom: circles show the measured radial velocities and 90% 
uncertainties. Stellar component radial velocities (solid and dashed lines in the lower panel) 



Herbig Sz Griffin 



and or bital phases were computed from the ephemeris and orbital solution of 
( I2OO6I ). Right panel: confidence countours (la, 2a, and 3a, inner to outer) for Doppler 



shift (Av/c) against the turbulent broadening parameter. These were computed for two 
groups of spectra, those nearest the minimum orbital radial velocity separation ( "low Av" ; 
dashed), and those nearest the maximum ("high Af"; solid). 
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Table 2. Line Flux Measurements 



OS Tmax Apred Ameas /meas fpred 

[logK] [A] [A] [10-6 phot cm^2 g-i] 



Fe XXV 


7 


8 


1 


8607 


1 


8684 


(0 


0073) 


3 


293 


(1 


846) 


1 


433 


Ca XIX 


7 


5 


3 


1772 


3 


1769 


(0 


0023) 


1 


135 


(0 


421) 





518 


Ca XIX 


7 


4 


3 


1909 


3 


1909 


(0 


0000) 





033 


(0 


223) 





172 


Ca XIX 


7 


5 


3 


2110 


3 


2110 


(0 


0000) 





373 


(0 


355) 





156 


Ar XVIII 


7 


7 


3 


7338 


3 


7355 


(0 


0081) 





328 


(0 


351) 





236 


Ar XVII 


7 


4 


3 


9491 


3 


9582 


(0 


0059) 





961 


(0 


558) 





397 


S XV 


7 


3 


4 


0883 


4 


0883 


(0 


0000) 





093 


(0 


226) 





049 


S XVI 


7 


6 


4 


7301 


4 


7272 


(0 


0026) 


1 


615 


(0 


530) 


1 


030 


Si XIV 


7 


4 


4 


9468 


4 


9468 


(0 


0000) 





193 


(0 


311) 





124 


S XV 


7 


2 


5 


0387 


5 


0405 


(0 


0035) 


1 


253 


(0 


556) 


1 


136 


S XV 


7 


2 


5 


0648 


5 


0648 


(0 


0000) 





541 


(0 


499) 





246 


S XV 


7 


2 


5 


1015 


5 


1038 


(0 


0051) 





991 


(0 


505) 





373 


Si XIV 


7 


4 


5 


2174 


5 


2172 


(0 


0065) 


1 


371 


(0 


604) 





580 


Si XIII 


7 


1 


5 


4045 


5 


4045 


(0 


0000) 





031 


(0 


255) 





131 


Si XIII 


7 


1 


5 


6805 


5 


6698 


(0 


0030) 


1 


507 


(0 


511) 





391 


Si XIII 


6 


9 


5 


8160 


5 


8042 


(0 


0035) 





815 


(0 


410) 





027 


Si XIV 


7 


4 


6 


1831 


6 


1825 


(0 


0008) 


4 


908 


(0 


466) 


4 


143 


Si XIII 


7 





6 


6480 


6 


6469 


(0 


0011) 


3 


341 


(0 


426) 


2 


937 


Mg XII 


7 


2 


7 


1063 


7 


1028 


(0 


0018) 


1 


095 


(0 


290) 





735 


Mg XI 


6 


9 


7 


3101 


7 


3101 


(0 


0000) 





173 


(0 


211) 





044 
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Table 2 — Continued 



Ion log Tmax Apred Ameas fmeas /pred 

[logK] [A] [A] [10-6 phot cm-^s-^] 



Fe XXII 


7 


1 


7.6812 


7.6812 


(0.0000) 





026 


(0 


136) 





073 


Al XII 


7 





7.7573 


7.7668 


(0.0054) 





521 


(0 


281) 





458 


Mg XI 


6 


9 


7.8503 


7.8503 


(0.0000) 





472 


(0 


277) 





257 


Al XII 


6 


9 


7.8721 


7.8721 


(0.0000) 





608 


(0 


300) 





319 


Fe XXIII 


7 


2 


7.9009 


7.9009 


(0.0000) 





079 


(0 


181) 





149 


Fe XXIV 


7 


4 


7.9857 


7.9794 


(0.0138) 





515 


(0 


484) 





523 


Fe XXIV 


7 


4 


7.9960 


7.9894 


(0.0121) 


1 


055 


(0 


733) 





266 


Fe XXIII 


7 


2 


8.3038 


8.3055 


(0.0025) 


1 


078 


(0 


346) 





552 


Fe XXIV 


7 


4 


8.3161 


8.3217 


(0.0022) 





901 


(0 


335) 





581 


Fe XXIV 


7 


4 


8.3761 


8.3810 


(0.0037) 





321 


(0 


272) 





225 


Mg XII 


7 


2 


8.4219 


8.4215 


(0.0008) 


5 


875 


(0 


531) 


5 


198 


Fe XXI 


7 


1 


8.5740 


8.5643 


(0.0056) 





377 


(0 


306) 





295 


Fe XXIII 


7 


2 


8.8149 


8.8158 


(0.0062) 





357 


(0 


302) 





566 


Fe XXII 


7 


1 


8.9748 


8.9741 


(0.0047) 





811 


(0 


360) 





582 


Mg XI 


6 


8 


9.1687 


9.1710 


(0.0015) 


2 


590 


(0 


475) 


1 


866 


Fe XXI 


7 


1 


9.1944 


9.1944 


(0.0000) 





249 


(0 


307) 





220 


Fe XXII 


7 


1 


9.3933 


9.3865 


(0.0150) 





277 


(0 


355) 





126 


Ne X 


7 





9.7083 


9.7082 


(0.0048) 


2 


002 


(0 


899) 


1 


210 


Ni XIX 


6 


8 


10.1100 


10.1169 


(0.0150) 





018 


(0 


220) 





044 


Fe XX 


7 





10.1203 


10.1332 (0.0032) 





654 


(0 


432) 





194 
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Table 2 — Continued 



OS Tjjiax Apred Ameas /meas /pred 

[logK] [A] [A] [10-6 phot cm-^s-i] 



Ne X 


7 





10 


2390 


10 


2390 


(0 


0015) 


4.084 


(0.634) 


3.736 


Fe XXIV 


7 


4 


10 


6190 


10 


6166 


(0 


0041) 


3.659 


(1.812) 


3.792 


Fe XIX 


6 


9 


10 


6491 


10 


6435 


(0 


0086) 


0.759 


(0.764) 


0.237 


Fe XXIV 


7 


4 


10 


6630 


10 


6650 


(0 


0027) 


2.314 


(0.652) 


1.990 


Fe XXIII 


7 


2 


11 


0190 


11 


0188 


(0 


0043) 


2.517 


(1.651) 


1.944 


Fe XXIV 


7 


4 


11 


0290 


11 


0312 


(0 


0039) 


3.934 


(1.563) 


2.460 


Fe XXIV 


7 


4 


11 


1760 


11 


1764 


(0 


0015) 


6.469 


(0.946) 


4.445 


Fe XVIII 


6 


8 


11 


3260 


11 


3132 


(0 


0041) 


1.337 


(0.574) 


0.469 


Ne IX 


6 


6 


11 


5440 


11 


5436 


(0 


0060) 


0.578 


(0.519) 


0.560 


Fe XXIII 


7 


2 


11 


7360 


11 


7400 


(0 


0012) 


7.656 


(1.068) 


6.279 


Fe XXII 


7 


1 


11 


7700 


11 


7719 


(0 


0014) 


6.244 


(0.996) 


5.384 


Ne X 


6 


9 


12 


1348 


12 


1352 


(0 


0010) 


30.490 


(2.300) 


25.433 


Fe XXIII 


7 


2 


12 


1610 


12 


1610 


(0 


0000) 


4.296 


(1.306) 


3.492 


Fe XVII 


6 


7 


12 


2660 


12 


2713 


(0 


0053) 


2.024 


(0.937) 


0.900 


Fe XX 


7 





13 


3850 


13 


3844 


(0 


0075) 


1.416 


(1.071) 


0.962 


Fe XIX 


6 


9 


13 


4230 


13 


4230 


(0 


0000) 


1.162 


(1.131) 


0.389 


Ne IX 


6 


6 


13 


4473 


13 


4493 


(0 


0050) 


5.083 


(2.097) 


3.874 


Fe XIX 


6 


9 


13 


4620 


13 


4620 


(0 


0000) 


0.893 


(1.591) 


0.883 


Fe XIX 


6 


9 


13 


4970 


13 


4970 


(0 


0000) 


2.554 


(1.311) 


1.552 


Fe XIX 


6 


9 


13 


5180 


13 


5324 


(0 


0048) 


4.074 


(1.938) 


3.422 
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Table 2 — Continued 



Ion log Tmax Apred Ameas /meas /pred 

[logK] [A] [A] [10-*^ phot s-i] 



Fe XIX 


6 


9 


13 


6450 


13 


6599 


(0 


0150) 


0.609 


(0.889) 





543 


Ne IX 


6 


6 


13 


6990 


13 


6927 


(0 


0043) 


3.873 


(1.503) 


1 


714 


Fe XIX 


6 


9 


13 


7950 


13 


7986 


(0 


0109) 


1.971 


(1.789) 


1 


361 


Fe XVII 


6 


7 


13 


8250 


13 


8351 


(0 


0052) 


3.759 


(1.585) 





718 


Fe XVIII 


6 


8 


14 


2080 


14 


2132 


(0 


0037) 


6.378 


(1.987) 


5 


219 


Fe XVIII 


6 


8 


14 


2560 


14 


2560 


(0 


0000) 


0.442 


(1.192) 





996 


Fe XX 


7 





14 


2670 


14 


2670 


(0 


0000) 


3.975 


(2.142) 


1 


300 


Fe XVIII 


6 


8 


14 


3430 


14 


3430 


(0 


0000) 


0.323 


(0.727) 





618 


Fe XVIII 


6 


8 


14 


3730 


14 


3730 


(0 


0000) 


1.205 


(1.127) 


1 


325 


Fe XVIII 


6 


8 


14 


5340 


14 


5355 


(0 


0050) 


3.860 


(1.546) 


1 


003 


Fe XVII 


6 


7 


15 


0140 


15 


0109 


(0 


0022) 


10.500 


(2.349) 


8 


749 


Fe XIX 


6 


9 


15 


0790 


15 


0825 


(0 


0068) 


2.224 


(1.393) 


1 


098 


VIII 


6 


7 


15 


1762 


15 


1680 


(0 


0075) 


1.598 


(1.276) 





664 


Fe XIX 


6 


9 


15 


1980 


15 


2082 


(0 


0058) 


2.238 


(1.398) 





928 


Fe XVII 


6 


7 


15 


2610 


15 


2600 


(0 


0061) 


2.795 


(1.488) 


2 


462 


Fe XVIII 


6 


8 


15 


8700 


15 


8572 


(0 


0039) 


2.800 


(1.661) 





432 


Fe XVIII 


6 


8 


16 


0710 


16 


0814 


(0 


0058) 


2.960 


(1.796) 


1 


813 


Fe XVII 


6 


7 


16 


7800 


16 


7650 


(0 


0029) 


3.114 


(2.137) 


3 


792 


Fe XVII 


6 


7 


17 


0510 


17 


0547 


(0 


0041) 


10.810 (3.391) 


4 


557 


Fe XVII 


6 


7 


17 


0960 


17 


1011 


(0 


0048) 


6.809 (2.930) 


4 


086 
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Table 2 — Continued 



Ion 


logTmax 


Apred 


X ^ 


f ^ 

J meas 


/pred 




[logK] 


[A] 


[A] 


[10~^ phot cm" 




Fe XVIII 


6.8 


17.6230 


17.6080 (0.0025) 


4.326 (2.586) 


1.323 


VII 


6.4 


17.7680 


17.7680 (0.0000) 


0.300 (0.315) 


0.041 


VII 


6.4 


18.6270 


18.6270 (0.0000) 


0.900 (0.926) 


0.123 


VIII 


6.7 


18.9698 


18.9744 (0.0043) 


17.680 (5.846) 


13.292 



Note. — These are the hues which were used in the emission measure and abun- 
dance reconstruction. The fpred, is what the emission measure and abundance model 
predicts. 

^Values in parentheses are la uncertainties. If the uncertainty for the wavelength 
is 0.0, then the line position was frozen in the fit. 
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We have applied this technique to several other spectra (e.g., see 



Huenemoerder et al 



20071 . and references therein). Given the form of the problem there is no unique solution. 
However, results can be useful for comparison of emission measures derived with similar 
methods. Figure [3] shows our reconstructed emission measure distribution whose values are 
also given in Table [3l The corresponding abundances are shown in Figure H] and are listed 
in Table H 
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Table 3. Emission Measure Model 



logT EM EMio^ EMhigh 

[log K] [lO^^cm-^] 



6.4 


2 45e-02 


4 60e-03 


1.31e-01 


6.5 


3 76e-02 


8 90e-03 


1 59e-01 


6 6 


4 40p-02 


1.41e-02 


1 37e-01 


6 7 


6 5Qp-n2 


3 00p-n2 


1.45e-01 


6.8 


1.74e-01 


1 03e-01 


2 96e-01 


6.9 


6.46e-01 


4.31e-01 


9.69e-01 


7.0 


1.53e+00 


l.lOe+00 


2.11e+00 


7.1 


1.09e+00 


7.33e-01 


1.63e+00 


7.2 


1.35e+00 


9.38e-01 


1.95e+00 


7.3 


2.05e+00 


1.29e+00 


3.26e+00 


7.4 


1.79C+00 


1.07C+00 


2.98e+00 


7.5 


l.lOe+00 


6.21e-01 


1.94e+00 


7.6 


7.48e-01 


3.35e-01 


1.67e+00 


7.7 


6.13e-01 


1.85e-01 


2.02e+00 


7.8 


4.65e-01 


l.OOe-01 


2.16e+00 


7.9 


2.47e-01 


4.79e-02 


1.27e+00 


8.0 


8.67e-02 


2.09e-02 


3.60e-01 
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Note. — The reconstructed emission mea- 
sure over the temperature range of sensitive 
features. The "low" and "high" values are the 
logarithmic la boundaries from Monte-Carlo 
iterations. 
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T ' ' ' ' 1 ' ' 1 ' 1 ' ' ' ^ 




Log 



Fig. 3. — Reconstructed emission measure, for a distance of 450 pc and Nh = 2 x lO^^cm"^. 
The upper and lower gray boundaries are la (logarithmic) statistical uncertainties from 
Monte-Carlo iterations in which the measured fluxes were perturbed randomly according to 
their measurement uncertainties. The integrated emission measure is 1.2 x lO^^cm"^. 



2.4- 1- Dynamics 

We searched for dynamical effects by using the mean plasma model spectrum as a 
template for fitting narrow regions by adjusting the Doppler shift, turbulent broadening 
velocity, and local normalizations. The 10-12.5 A region has a number of lines which 
make it useful for this purpose, and we used that entire wavelength interval in the fits. 
We fit each observation independently for line-of-sight velocity. The results are listed 
in Table [1] and shown in Figure O The radial velocities are consistent with the orbital 
velocities of the stellar components — that is, they are generally less than or equal to the 
orbital radial velocities, but are limited by the spectrometer's sensitivity. According to the 
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Table 4. Relative Elemental Abundances 



Element Abund^ FIP'^ 

[eV] 






0.06 


Ne 


0.42 


Mg 


0.19 


Al 


0.64 


Si 


0.17 


S 


0.12 


Ar 


0.20 


Ca 


0.55 


Fe 


0.13 


Ni 


0.14 



0.03) 


13.618 


0.03) 


21.564 


0.02) 


7.646 


0.28) 


5.986 


0.01) 


8.151 


0.03) 


10.360 


0.17) 


15.759 


0.33) 


6.113 


0.01) 


7.870 


0.14) 


7.635 



Note. — Uncertainties in "()" are Icr values. 



Ab undances are relativ e to So lar using the val 



ues of 



Anders fc Grevessd (119891 ) 



TIP is the first ionization potential. 



-24- 




Fig. 4. — 9^ Ori E abundances (circles) relative to Solar photospheric values 



flAnders &: Grevesse 



(Cunha. Hubenv & Lanz 


Cunha & Smith 


2005 


) or 



19891 ) as derived from emission measure reconstruction. Error-bars 

give Orion stellar photospheric 



2006 



Cunha &: Lambert 



1994 



Cunha. Smith &: Lambert 



1998 



20041 ) abundance ratios for Fe, Si, S, O, 



Ar, and Ne. 



orbital solution of iHerbig &: GrifEnI (120061 ). the stellar components have a projected radial 
velocity amplitude of about 80kms~^ and a systemic velocity of 30kms~^ (note that for 
the adopted masses of S.SM©, the inclination is i = 61°). While there is a slight systematic 
offset towards the red-shifted component at some phases (e.g., 0.28, 0.30, 0.33, 0.35), the 
trend does not persist since at other phases (e.g., 0.55, 0.60), measurements span the range 
between the components' orbital velocities. 

To obtain better sensitivity and to examine line widths, we also fit spectra combined 
into two groups, one nearest to zero velocity separation (near phase 0.25; 4 observations) and 
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the other nearest to maximum velocity separation (near phase 0.5 or 1.0; 6 observations). 
We computed contour maps in hne-of-sight and turbulent broadening velocities. The 
point of this is not that we necessarily expect turbulent broadening, but that there could 
be broadening due to binary orbital effects, and fitting turbulent broadening is simply a 
useful parameterization of this. For instance, for equally X-ray bright stars, the X-ray 
line's measured radial velocity could always be zero (or the systemic value), but the lines 
could broaden and narrow, modulated by the orbital radial velocity. If only one star 
were the X-ray source, the lines could shift but maintain constant width. Since we are 
photon-limited, we need to group spectra in order to obtain significant counting statistics. 

The resulting confidence contours are shown in Figure [2l We see the radial velocity 
offset towards a small positive velocity as we should in the low-Af group at a level of 
about 30-120 kms~^ (or Af/c ~ 1-4 x 10~^; 90% confidence limits), while the high-Af 
group range is —60 to 60kms~^. Broadening is marginally significant with 90% confidence 
contours of 60 - 300 km s^^ for the high-Af group and - 200kms~^ for the low-Af group 
(with a la lower limit of 20kms~^). 



2.4-2. He-like Triplets 



The HETGS bandpass includes the He-like triplet lines of Mg xi, Ne ix, and O Vll, 



1969 



19731). 



which are useful diagnostics of density in the coronal regime f lGabriel &: Jordan 
Due to the absorption towards Orion and the low sensitivity of the HETGS at 22 A, we 
have no useful data on O vil, but we do have spectra of Mg xi and Ne ix. By fitting the 
line ratios of the resonance (r), intercombination (z), and forbidden (/) lines we are able to 
put upper limits on the coronal electron density. In our fits, we constrained the relative 
positions of the triplet components and included blends as estimated from the emission 
measure model spectrum. The Ne x Lyman series converges near the Mg xi i-line (9.230 
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A). We included relevant lines of this series by providing an initial guess for their strengths 
by scaling fluxes according to their relative /-values from the isolated and well detected 
Ly-7 and 6 lines, since APED does not contain lines with upper levels n > 5. The locations 
of the weaker and blended features are 9.215 A (Ne x Ly-a series n = 9 to n = 1 transition), 
9.246 A (Ne x 8-to-l transition), and 9.194 A (Fe xxi). The weakest (Fe xxi) line's position 
was frozen, while the others positions and fluxes were left free and gave reasonable fitted 
wavelengths. While this is not a complete and accurate plasma model for the region, it is a 
useful parameterization to obtain ratios for the interesting features. 

There are Fe blends in the Ne ix triplet region, but since the Ne:Fe abundance ratio is 
high (see Table H]), these are not severe. 

The continuum in each region was evaluated from the plasma model and not governed 
by free parameters. 

Results for both triplets are consistent with low density. In Figure [5] we show the 
spectra and confidence contours of the G and R ratios, defined as G = {f + i)/r and 



R = f /i. G is primarily a function of temperature, a nd R of density. T he density and 

)£1 



temperature dependence are from APEC calculations (iSmith et al. 



20011) 



^Data for the lines are available from|http : //cxc . harvard . edu/atomdb/f eatures/denHETG . ps 
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Fig. 5. — The upper plots show the Mg xi (left) and Ne ix (right) triplets counts spectra 
(black histogram), a multi-Gaussian fit (gray line) and residuals (lower panels of top plots). 
The positions of the triplet components are marked. The lower plots show the fits to the G 
(= (/ + i)/r) and R {= f /i) ratios, which are primarily temperature and density sensitive, 
respectively, as given by the grid. The contours are the Icr, 2cr, and 3(7 levels from lower /inner 
to upper/outer contours. 
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2.4-3. Abundance Ratios 

By forming appropriately weighted flux ratios of line pairs, we can obtain relatively 
temperature-insensitive abundance ratios. This is achieved by relating a linear combination 
of emissivities of He and H-like resonance lines for one element to a similar linear 
combination of another element, so as to minimize fluctuations with temperature in 
their ratio. If this ratio is constant, then the analogous combin ation of line f luxes gives 



the abundance ratio, 
described generally by 



his technique was explaine d in detail by iLiefke et al. 



Garcia- Alvarez et al 



fl2005h . 



2008|) and 



We define the ratio in the following way: 

_ + aiFi,2 _ Aiao / [ei,i(T) + aiei,2{T)]D{T) dT 



(1) 



i"2,i + 02^2,2 ^2 /[e2,i(T) + a2e2,2(T)]D(T)rfT 
in which subscripts i,j on F refer to the line j from element i, F is the observed flux, e(T) 
is the emissivity, D{T) the emission measure, and Ai the abundance of element i. The 
parameters, a„ are to-be-determined. It is clear that if the terms in square brackets within 
the integrals are identical in numerator and denominator for all T, then the integrals cancel 
and we are left with the abundance ratio (parameters Oi and 02 serve to flatten the ratio, 
and Oo normalizes it). We determined the parameters by minimizing the variation in the 
ratios using H-like and He-like line emissivities from the APED database. Re-writing in 
terms of the abundance ratio for these lines, we have 

^2 ~ Uo/ F2,H + a2F2,He 

in which subscripts H and He represent the hydrogen-like Lyman-a doublet and the 
helium-like Lyman-a resonance line, respectively. We tabulate coefficients for a few useful 
ratios in Table [5] as derived from APED for units of F in [phot cm^^s^^] and abundances 
relative to Solar. Minimization of the ratio variances were restricted to temperature ranges 
where the emissivities were greater than 1% of their maximum. In Table [6] we give the 
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abundance ratios derived from the temperature-insensitive ratios along with values from 
emission measure modeling. The ratios from each method are in very good agreement. This 
means that abundance ratios can be derived fairly easily, without resorting to emission 
measure reconstruction. 
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Table 5. Temperature Insensitive Ratio Coefficients 



Ratio 






^2 


a 


Ne:Mg 


1.440 


0.050 


2.675 


0.273 


Mg:Si 


0.536 


0.085 


1.920 


0.064 


Si:S 


1.197 


0.125 


1.545 


0.092 


S:0 


0.212 


2.990 


0.000 


0.137 


Ne:0 


0.347 


2.590 


0.000 


0.120 



Note. — The coefficients are to be used as defined by 
Equation[2l The last column, a, is the standard deviation 
of the emissivity ratio, which is never perfectly fiat, and 
coefficient Oq represents the mean of the ratio. 
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Table 6. Abundance Ratios' 



Elements 


^ (Tinsens) 


r{EM) 


Ne:Mg 


1.7 (0.3) 


2.2 (0.2) 


Mg:bl 


i.U [V.i ) 


1.1 (0.1 ) 


Si:S 


1.2 (0.2) 


1.4 (0.4) 


S:0 


1.4 (1.0) 


1.9 (0.8) 


Ne:0 


4.9 (2.8) 


6.5 (2.3) 


Mg:0^ 


2.9 (1.8) 


2.9 (1.0) 


Si:0^ 


1.8 (0.9) 


2.6 (0.9) 


Note. — 




gives the abundance ratios from 



the temperature-insensitive method, while r{EM) gives 
ratios from the emission measure reconstruction (see Ta- 
ble Hj). Values in parentheses are la statistical uncertain- 
ties. 

^The ratios are of relative Solar photospheric abun- 
dances. 

^Derived from preceding ratios - not determined di- 
rectly. 
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3. Discussion, Interpretation 



The recent determination by lHerbig fc GriffinI (120061 ) that 6^ Ori E is a moderate mass 
pre-main-sequence spectroscopic binary is very important in the context of stellar evolution 
and X-ray activity. When 6^ Ori E arrives on the main sequence, we expect it to be faint 
or non-detectable in soft X-rays. Yet at the age of 0.5 Myr, it is the second-brightest steady 
X-ray source in the Orion Trapezium . The binary system ha s = 1.2 x 10^^ ergs s~^, and 



given an optical luminosity of 29L0 (IHerbig fc Griffin 



2006f ) it thus has L^/L 



for t 



10" 



le pair. This value is near the saturation limit of coronally active stars ( iProsser et al. 
19961 ). The X-ray emission is similar to other magnetically active stars, having a broad 
temperature distribution and narrow emission lines. Hence we surmise that 6^ Ori E has 
dynamo activity and probably has stron g convection zones. This is also in accordance with 



evolutionary models of stellar interiors (ISiess. Dufour fc Forestini 



2000l ) which indicate a 



substantial convec tion zone for stars like the 6^ Ori E components (also see Figure 1 of 



Stelzer et al. 



20051). 



Prior to the 



Herbig fc GriffinI (j2006|) deter mination that 6^ Ori E is a binary of G-type 



stars, 6^ Ori E was considered to be a B5-star. 



Schulz et al. 



(120031) considered the X-rays 



to be a hybrid of wind shock emission and n iagnetically c o nfined winds, but they did note a 



striking similarity to active coronal sources. 



Stelzer et al. 



(120051 ) interpreted the emission as 



from a weak wind, but noted unusual emission characteristics, such as due to an extended 
magnetosphere and magnetically confined wind shocks. It is now clear in hindsight, given 
the spectral types, that emission is coronal in nature. 

6^ Ori E did not show any distinct fiares during our observations, which is somewhat 
unusual for coronally active stars. Th is lack of activity is consistent, however, with the 



long intervals of constant fiux seen by 



Stelzer et al. 



(120051 ). There was one observation in 



which the fiux was /ow;er than our average (see Figure [2]), and examination of the spectra 



- 33 - 



shows that this is manifested in diminished short wavelength flux (below about lOA). This 
variation between observations cannot be attributed to rotational phase dependence since it 
does not repeat. We also found no significant variability within any of our observations. It 
is possible that the system is so coronally active that a significant proportion of the average 
flux is from continuously visible flares, which would also give rise to the dominant emission 
measure peak at logT = 7.3, but one observation had a bit less flaring. 



Stelzer et al. 



(I2OO5I ) from heavily 



The HETGS flux was about half that reported by 
piled, low resolution spectra. The HETGS flux calibration is accurate to about 5%. 
Since flux for our observations was effectively constant and spans the time of the COUP 



observations, the difference in flux without obvious flares is unusual. Since the 



Stelzer et al. 



( I2OO5I ) analysis was made difficult by the high photon pile-up in the core of 6^ Ori E, they 
resorted to spectral analysis of photons only from the wings of the point-spread-function. 
We reanalyzed spectra for 6^ Ori E from one of the COUP datasets, observation ID 4373, 
for which the flux was constant. We used an extraction radius of 2.25 arcsec centered on 
the source, including piled photons and made standard responses for the region. To fit the 



spectrum, we used the pileup model of 



David (I2OOII ) as implemented in ISIS. Pileup is a 



very non-linear process; there can be multiple solutions since the count-rate first saturates 
with increasing fluence, then can decrease as events are rejected from telemetry. Finally, for 
extremely high pileup, the count rate can again increase when the core is fully saturated 
and the wings grow. We used two-temperatu re component, absorbed APED plasma 
models, similar to those of 1 



Stelzer et al. 



( 120051 ) ■ and used Monte-Carlo techniques to explore 



parameter space, given the multi-valued nature of pileup fitting and the possibly degenerate 
nature of the models. While we found solutions with fluxes similar to those presented by 



Stelzer et al 



(l2005l ). we found equally acceptable solutions (reduced < 1-3 ) with fluxes 
comparable to our HETGS-derived values. Incidentally, all our fits to this one spectrum 
preferred > 3 x lO^^cm"^, a bit larger than our adopted value of 2 x lO^^cm"^. We 



-34- 



conclude that the IStelzer et aLl (120051 ) flux is probably in error and the source is probably 
steady outside of distinct flares with a flux of about 3 x 10^^^ ergs cm~^s~^. 

The radial velocities determined from lines are consistent with the orbital dynamics. 
Given a peak-to-peak orbital radial velocity amplitude of 160kms~^, we have marginal 
sensitivity for detection of the orbital modulation if emission were dominated by one stellar 
component (see the confidence limits in Tableland Figure [2]). Variability over the time 
period of observations could destroy any orbital systematic radial velocities in X-ray lines 
if the relative activity level of the two stellar components changed. Our phase coverage is 
also poor, but the lack of significant velocity offsets at phases of maximum orbital velocity 
separation suggests that both stellar components are roughly equal in X-ray emission. The 
marginal detection of line broadening, particularly at these same phases, is consistent with 
the broadening being due to orbital velocity effects. We conclude that the lines are similar 
to other coronal sources - narrow, and effectively unresolved. 

The absolute abundances are rather low when compared to other coronal sources. 
If we compare Ori E to the abundances derived from low-resolution COUP spectra 



of 



Maggio et al. 



(120071 ). they are not only lower by a factor of 5 or more in general, the 
ratios are also different — they are, in fact, uncorrelated. This probably has as much to 
do with different methods and spe ctral resolutio n s tha n with intrinsic differences between 



9^ Ori E and average Orion stars. 



Maggio et al 



(120071 ) use two temperature component 



fits, and these cannot accurately reproduce abundances and emission measures for realistic, 
continuous emission measure distribution plasmas. If a fitted temperature component is off 
the peak of some ion's temperature of peak emissivity, and there is actually plasma at that 
temperature, then a two-temperature model will artificially increase the abundance of that 
element in order to reproduce the flux. Comparison of low and high resolution results with 
different modeling approach is in general not meaningful. 
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3.1. Loop Sizes 

The geometric structure of stellar coronae is largely an open question, and is relevant 
to energetics, variability, and likelihood of interactions with stellar companions or disks. 
While there are many uncertain parameters, we can provide order-of-magnitude estimates 
via several methods to show that loops could be compact (small fraction of the stellar radii 
of 7R0; see Table [7]), or comparable to the stellar radii an d thus a significan t fract ion of the 



stellar separation (the semi-major axis is about 2.5-R*; see IHerbig fc Griffiru (1200 



3)). 



If we assume that the X-ray emission originates in an ensemble of identical semi-circular 
loops, we can estimate the order of magnitude of the loops' radius. The loop radius (or 
height if vertically oriented) relative to the stellar radius can be expressed as 

/z = E,f i?-n-/^iV-/Vf (3) 

in which E is the volume emission measure, R the stellar radius, N the number of identical 
loops, n the electron density, and a is the loop aspect ratio (cross sectional radius to 
height, < 1), and the subscripts indicate the power of 10 scale factor (all cgs or unit-less 
quantities). Only two of these parameters are well determined, E = 1.2 x 10^^ cm~^ from 
ou r X-ray spectral r nodeli ng, and i? ~ 5 x 10^^ cm from the radial velocity curve analysis 



of iHerbig fc GrifEnI (120061 ). Densities are poorly constrained; we will adopt lO^^cm"^ for 
argument (see Figure [5]). From Solar loops, a is about 0.1, and given the lack of variability 
in 6^ Ori E (to about 10% accuracy; see Figure [21), we will let = 100. If we assume that 
the emission is divided equally between the two binary stellar components, then with these 
parameters we obtain h ~ 0.02, implying that the coronae are compact. For a density of 
10^°cm~^, this height increases by a factor of 20 to about 0.4, a significant fraction of the 
orbital separation. 

We can also estimate loop parameters using hydrodynamical models from the flare 
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temporal and spectra l properties. Assuming that a single flaring loop dominates the 
emission 



Reald (l2007l ) expressed the loop size (his Equation 12) as Lg ~ 3(To/T/vf)^ Tq(^ tm, 



3- 



Here Lg is the loop half-length in units of 10^ cm, Tq is the maximum temperature during 
the flare (and Tqj is the same in units of 10^ K), Tm is the temperature at which maximum 
density occurs, and tm,3 is the time from flare start (in ks) at which maximum density 
occurs. To apply such a model in detail, we would need the evolution of emission measure 
(a proxy for density) and temperature (from time-resolved spectra) from the rise to the 
decay of a flare. We do not have such, so we will make some reasonable approximations to 
obtain an order-of-magnitude hydrodynamical loop size. 



2005r i: from this. 



The COUP observation of 6^ Ori E detected a flare (jStelzer et al.l 
we estimate that Tm,3 = 25. We have an emission measure distribution; we assume that 
the hotter peak and hot tail represents the integrated history over many flares. From 
time-resolved analyse s of other stars, we have seen that such a hot peak can be directly 



attributed to flaring (IHuenemoerder. Canizares fc Schulz 



2001 



Giidel et al. 



20041 ). We will 



thus assume that the flare mean temperature of maximum density (or maximum emission 
measure) corresponds to our strongest peak, or log T^f = 7.3 (see Figure [3]). We will assume 
that the hot tail of the EMD represents the maximum flare temperature. This is less well 
deflned, and we adopt logTo = 7.8, which is where there is an inflection in our EMD. From 
these parameters we flnd a relative loop half-length of about 3.8 stellar radii, or a height 
(for semi-circular vertical loops) of 2.4 stellar radii. 



We can also adopt flare parameters from other giant stars, such as HR 9024 (ITesta et al. 



2003) (also see Table [ID, for which logTo = 7.9 and Tq/Tm = 1.4. Thus, if 6^ Ori E flare 
temperatures and densities are similar to those on HR 9024, the loop half-length is about 1 
stellar radius (or a height of 0.6 stellar radii). 



Since these are order-of-magnitude estimates, we conclude that flare loops can be of 
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order the stellar radius. In sum, an origin of the emission from magnetically confined coronal 
loops is not unreasonable. To better determine the coronal geometry, more information is 
needed, such as more stringent constraints on density, detection of rotational modulation, 
or time-resolved spectroscopy of a large flare. 



4. Comparison to Other Stars 



To understand the nature of the X-ray emission from 6^ Ori E, we must examine it 
in the evolutionary context of stars of similar mass. We have collected information for 
several other stars, both pre-main-sequence and post-main-sequence, with masses ranging 
from 2 to 3.5M0. Information and sources are listed in Tableland in Figure [6] we show 
the objects on a temperature-luminosity diagram along with evolutionary tracks. Of the 
sample, 6'^ Ori E has the highest relativ e X-ray luminosity, being as high in L^/Lhni as 



saturated" short p eriod active binaries (IVilhu fc Rucinski 



Prosser et al. 



1983 



Cruddace fc Dupree 



1984 



19961 ). even though its period is somewhat longer than those systems. If we 



look to the future of 6^ Ori E's evolution and consider AB Aur, we see that L^/ L]^o\ may 
become much smaller; we expect main-sequence A-type stars to be very faint or undetected 
in X-rays (for reference, a spectral type AO V star has a mass of ~ 2Mq, and a B5 V of 
~ 6M0). 

Stellar rotation is well known to be a key factor in magnetic dynamo generation. If 
we compare the sample's X-ray activity as a function of period (Figure [Tj), excluding AB 
Aur, we see a strong anti-correlation which holds for both pre- and post-main-sequence 



objects. This is similar to the behavior o: 



stellar coronae OWalter fc Bowyer 



1981 



active giants, binari e s, or main sequ ence late-type 



Pizzolato et all 12003 



Gondoin 



2005h . AB Aur, 



in spite of its short period, has very low activity; it is approaching the low-activity 
main-sequence era of its life, and may have a radically different emission mechanism, such 
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as fro m wind or accretion affects, with lower Lx and X-ray temperatures (iTelleschi et al 



20071 ). It's convective zone, necessary for magnetic dynamo generation, is quite small. 



being less than 1% of the stellar radius, compared to about 20% for the other stars 



(ISiess. Dufour fc Forestini 



2000 



Stelzer et al. 



2005 



their Figure 1), so it is reasonable to 



exclude it from period-activity relations of stars with significant convective regions. 

9^ Ori E has a very hot corona, characterized by an emission measure distribution 
with a strong peak at about 20 MK (see Figure [3]). It is similar to SU Aur and 
HR 9024 (see Table [7]), two objects which displayed strong flares during their X-ray 
observations. A hot emission measure peak has been directly identified with flares 



(IHuenemoerder. Canizares fc Schulz 



2001 



Giidel et al. 



20041 ). In this context, it is 



curious that during all the COUP and HETG exposures, only one distinct flare was seen 



( Stelzer et al. 



20051 ) . We can only speculate that perhaps 9^ Ori E is so active that the flare 



rate is so high that they are nearly always superimpos ed and create a nearly co nstant flux 



Caramazza et al 



(120071). Our single 



Such was found plausible for Orion's low-mass stars by 
low-flux HETGS observation (see Figure [2]) did have a relative deficit in short wavelength 
flux (< 5 A), a region sensitive to the highest temperature plasmas; it is consistent with 
diminished flare activity. To sustain continuous flaring, there has to be a continuous source 
of erupting magnetic fields and their reconnection. The pressure scale-height for a hot, 
low-density plasma is a significant fraction of the binary stellar separation of 9^ Ori E. 
There could be star-to-star magnetic reconnections sustained at a fairly high level by the 
binary proximity and dynamo action generating sufficiently large loops. 

The presence of a very hot corona and the low probability of distinct flares is com mon 



to several evolved giants, such as HR 9024, u Vel, 31 Com, or 



2007 



Testa. Drake fc Peres 



2004 



Ayres. Hodges-Kluck fc Brown 



M Peg (ITesta et al. 



20071 ). While this is not 



understood, it could be a significant trait of the coronal heating mechanisms. 
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Another distinguishing characteristic of 6^ Ori E is the relatively low mean metal 
abundance. Table H] shows relative elemental abundances from the EMD analysis. All 
abundances are significantly below unity with oxygen at an extrem ely low value. Wh en 



compared to abundances deduced in a similar analysis of 9^ Ori C (iSchulz et al. 



20031 ). its 



massive neighbor within the Orion Trapezium, then there are a few remarkable differences 
to note. Values for Ne, Al, and Ca seem very similar within uncertainties, while Mg, Si, 
S, and Ar are very different, being near or above unity in 6^ Ori C; values for O and 
Fe are even lower than in 6^ Ori C. Figure H] also compares coronal with average Orion 
stellar photospheric and nebular values. Since the photospheric and nebular abundances 
are all near unity, it seems that abundances from the hot X-ray plasmas are fundamentally 
different. A two-temperature modeling of Ori C could reconcile deficient O and Ne 



values by requiring a significantly higher column densi tv (IGagne et a. 



forcing the Ne/0 ratio to be similar to other coronae (iDrake fc Testa 



20051) . as would 



20051 ). though Fe 



would still remain low. In the case of 6^ Ori E, the application of a higher column will 
not change the low values for Mg and higher-Z elements, but also not enough for O, which 
at A{0) = 0.06 is extremely low. We also find that multi-temperature plasmas more 
accurately describe the spectra in both stars and thus see trends in Orion Trapezium stars 
which include some neon deficiency but clearly hot plasmas with low iron and oxygen 
abundances. For pre-main-sequence stars, depletion of metals has been explained by 
formation of dust in the disk, and the remainin g gas, seen heated to X-ray temp eratures in 



an accretion shock, is the metal-poor material (IDrake. Testa fc Hartmann 



20051 ). Magnetic 



coronae generally show low metals , but t end to have h igher Ne and O (e .g., II Peg 



Huenemoerder. Canizares fc Schulzl (120011 ). or HR 1099 



Drake et al. 



(1200 If ) I Hot plasmas 



of Orion's stars, if the two studied are representative, seem to be different. There is no 
theoretical explanation for coronal abundances, but these differences may be clues to 
fractionation mechanisms. 
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Table 7. Stellar Comparison 



Quantity 


0^ Ori E ^ 


SU Aur 


AB Aur 


// Vel 


HR 9024b CapellaAa^ 


Spec. Type. 


G8 III+G8 III 


G2 IV 


AO 


G5 III 


Gl III 


G8 III 


Teff [K] 


5012 


5550 


9750 


5030 


5530 


5000 


Age [My] 


0.5 


4 


4 


360 


post-ms 


525 


M/Mq 


3.5 


2.0 


2.7 


3 


2.9 


2.7 


R/Rq 


7 


2.6 


2.3 


14.4 


13.6 


12.2 


g/ge 


0.07 


0.3 


0.5 


0.01 


0.02 


0.02 


P [day] 


9.9 


1.8 


1.4: 


117 


23.2 


104 


vsini [kms~^] 


37 


66 


80 


6.2 


22 


5 


-^bol/L© 


29 


6.3 


49 


108 


90 


78.5 


L,/10=^0[ergss-i] 


121 


8 


0.4 


2.2 


63,125 


1.6 


(LjLboi)/10-6 


1100 


300 


2 


5.4 


181,360 


5.3 


VEM/IQ^"^ [cm-3] 


1200 


50 


5 


8 


200,800 


4 


[MK]'i 


20 


20 


4.7 


7.9 


30,80 


6.3 


A{FeY 


0.1 


0.6 


0.3 


1.8 


0.2,0.7 


0.8 


A{OY 


0.1 


0.3 


0.2 


0.7 


0.6,1.1 


0.4 


A{NeY 


0.4 


1.3 


0.6 


1.4 


1.1,1.2 


0.7 



^X-ray properties are for the binary system; if each stellar component contributes equally, 
values should be divided by two. 
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^When two quantities are given they are for quiescent and flare states, respectively. 
*^X-ray properties assume that component Aa dominates. 

'^A qualitative temperature of the high energy emission measure distribution, adopted 
from visual inspection of published curves or few-T fits. 

''Abundances are relative to Solar. 



Note. — Sources: 9^ Ori E: 


Herbie & Griffin 


( 


2006 


; this 




)aper. SU Aur: 


Franciosini et a 


I. (2007 


): Robrade & Schmitt ( 




DeWarf et al. 


( 


2003). AB Aur: 


Telleschi et al. 


(200"; 


). u Vel: I 


Wres, Hodees-Kluck & Brown 


(,2007i): (Wood et al. 


(2005). HR 9024: 


Avres. Hodees-Kl 


uck & Brown (2007 


Testa et al 


(2007). Capella: 


Ishibashi et al. 


(2006 


[ 


Canizares et a 


L (200ol): 


Hummel et al 


(199^ 


1); 


Ness et al. (2008): 


Gu et al. (2006^ 


. 
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o 
o 
o 



o 
o 




Te« [K] 



Fig. 6. — The stars listed in Table [7] are shown on a Temperature-Luminosity diagram. Pre- 
main-sequence stars are shown with open symbols and post-main-sequence stars are shown as 
solid symbols. The shapes encode L^jL^ai as a "star" for > 10~^, circle for > 10^'' to < 10^^, 



I Siess. Dufou 


r & Forestini 


Schaller et al. 


( 


1992 


). The 



(120001 ). The post-main-sequence tracks (thin gray lines) are from 
tracks are shown for 2, 3, and 4 . The post-main-sequence tracks were scaled slightly 
(0-3% in temperature, 8-12% in luminosity) to better coincide with the pre-main-sequence 
tracks.) 
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o 
o 
o 



o 



to 

o 



T 1 1 — I I I I 



o 

SU Aur 



o 

AB Aur 
J I 



T 1 1 1 — I I I I 



o 

ei Ori E 



o 

HR 9024-F 

O 

HR 9024-Q 



10 

P [day] 



H Vel 

CO 
Capella Aa 



_i I ' ' ' ■ 



100 



Fig. 7. — The normalized X-ray luminosity for stars listed in Table [7] are shown as a function 
of their periods. The stars largely follow an expected rotation-activity trend, except for 
AB Aur, which is an A-star near the main-sequence and is expected to be faint in X-rays. 
Lx/Lboi for 6*^ Ori E is for the system; if each stellar component contributes equally, its 
value should be divided by two. 
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5. Conclusions 

6^ Ori E is perliaps tlie only case known of a pre- main- sequence 3 Solar mass G-star 
binary. As such, it holds an important place in our understanding of X-ray dynamo 
generation and evolution. We believe that when it reaches the main sequence, it will emit 
a negligible fraction of its luminosity in X-rays. Yet now it is the second brightest X-ray 
source in the Trapezium. Furthermore, its relative X-ray luminosity (L^/Lboi) makes it as 
strong as any of the coronally active binaries. Thus we conclude that as moderate mass stars 
collapse toward the main sequence, they go through a phase of strong magnetic dynamo 
generation, very similar or identical to that of coronally active late-type stars which sustain 
a convective zone and shear-generated magnetic dynamo. Since A-stars are dark or at most 



quite faint in X-rays (ISchroder fc Schmitt 



20071 ). at some point, the dynamo vanishes, and 



the magnetic fields dissipate. AB Aur, which is near the main sequence and relatively faint 
in X-rays, is a possible future state of 6^ Ori E. It is also clear from the post-main-sequence 
objects that a dynamo can be generated as stars of this mass evolve into giants. 

Hot plasma abundances of 6^ Ori E (as well as of 6^ Ori C) are different from coronally 
active stars. 6^ Ori E is hot and has fairly steady X-ray emission. The high temperature is 
either due to unresolved flares, or to an unknown mechanism which also may be common 
to other active G-giants. The temperature structure, abundances, and low variability may 
be clues to plasma heating mechanisms. Continuing high-resolution spectroscopic studies of 
Orion stars will show us if some of these patterns are common in the Orion Nebular Cluster. 



Facilities: CXO (HETGS) 
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